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Q\ ' Why is non-Gaussianity interesting? One of generic predictions from inflation- 

ary scenarios is that primordial fluctuations are exactly Gaussian in linear or- 
der; however, the non-linearity in the inflation will produce weak non-Gaussianity. 
Thus, measuring the non-Gaussianity in the cosmic microwave background radi- 
' ation anisotropy is a probe of the non-linear physics in the very early universe. 

' Since the angular three point function is zero for the Gaussian field, it is sensitive 

I i to the non-Gaussianity. We predict its harmonic transform counterpart, the an- 

, gular Bispectrum, down to arcminutes angular scales, including the full effect of 

the radiation transfer function. We find that even the Planck experiment cannot 



o 
o 



detect the primary bispectrum from the inflation, as long as the single field slow- 
roll inflation is right. Non-linearities in the low redshift universe also produce the 
non-Gaussianity. We find that secondary bispectra are detectable by both MAP 
and Planck experiments. The secondary bispectra probe the non-linear physics of 
the low-redshift universe. Although this could be a contaminant to the primary 
signal, MAP and Planck experiments are found to be able to separate the primary 
^ ' from secondary effects well. We present a tentative comparison of the primary 

' bispectrum to the published COBE 4 year bispectrum. The data put a weak con- 

jy^ i straint on the parameter, and the constraint would become much tighter when we 

, use all modes available in the COBE data, and certainly forthcoming satellite ex- 

periments. As a conclusion, the bispectrum is a key measure to confirm or destroy 
the simple inflationary scenario in non-linear order that seems quite successful in 
" linear order. 

h . 

. F~ , 1 The CMB non-Gaussianity from inflation 



1.1 



Order- of -magnitude estimation of the non-linear coupling parameter /jvx 



We briefly consider the mechanism producing the primary CMB non-Gaussianity 
from the inflation. Since the observed CMB anisotropy AT/T is thought to be 
generated by curvature perturbations 4>, the primary non-Gaussianity should be 
encapsulated in <i>. The inflation makes weakly non-Gaussian. However, even if 
$ is Gaussian, AT/T can be non-Gaussian. 

AT/T is related to $ through the non-linear relation: 

AT 

— ~ g T ($ + / $ <E- 2 ) , (1) 

where gx is the radiatipn transfer function converting $ to AT/T. gr — —1/3 for 
the Sachs- Wolfe cffcctEl. The second term with /$ ~ 0(1) is the higher order cor- 
rection that arises from the second order gravitational perturbation theorya. Thus, 
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it makes AT weakly non-Gaussian, even when <3> is Gaussian. <I> is actually a weakly 
non-Gaussian field, and thought to be produced from the inflaton fluctuations 5<fi 
through the relational: 



$ ~ TO 



l g 9 [H + m-lf H 5(j> 2 y (2) 



g<s> ~ O(10), and fs$ ~ O(10 _1 ) for a class of slow-roll inflationary scenarios 
with the single inflaton field. If the dynamics of 8<j> field is simple enough to keep 
itself Gaussian throughout the evolution, then we can stop our consideration here. 
However, in generic, the non-trivial form of the interaction term in the equation of 
motiorn, or the non-linear coupling between the coarse-grained classical fluctuations 



and the fine scale quantum fluctuation^ make 84> weakly non-Gaussian, even if 
it is originally Gaussian field. The result is simply written as 

5^ ~ 9S4> [V + m~ l 1 f v rj 2 ^j , (3) 

where r\ is Gaussian, g$^ ~ 1, and f n ~ O(10~ 1 ). 
It is then useful to rewritten equations (Q)-^) as 



AT 

9t 



*l + (/* + .9$ 1 fs* + 9^94 fv) *l] . (4) 



T 

where $l = g^gs^t] ~ 10r/ is the Gaussian part of curvature perturbations. More 
specifically, we define the non-linear coupling parameter, Jnl, as 

$(x) = 4> L (x) + f NL ($|(x) - ($|(x)» , (5) 

where the angular bracket denotes the statistical ensemble average. Note that 
Jnl ~ /# + g^, 1 fstjj + g$> gjlfri> an d the first term ~ 0(1) from the second order 
gravity effect is dominant compared to other two, terms from slow-roll inflation 
~ O(10 -2 ). Jnl corresponds to a$ in Verde et al.0. 

1.2 The reduced angular bispectrum 

We use the angular bispectrum, harmonic transform of the angular three point 
function, as a measure of the CMB non-Gaussianity. 

AT(h)/T in the sky is expanded into the spherical harmonics: AT(n)/T = 
ozm^m(n). ai m is then written in terms of $-field, 

^)g Tl (k)Y; m (k). (6) 

$(k) is weakly non-Gaussian field given by the Fourier transform of equation (||). 
Instead of assuming the Sachs- Wolfe's —1/3 law for gri(k), we use the full gn{k) 
down to small scales by solving the Boltzmann transport equation with the CMB- 
FAST codela. We then define the reduced angular bispectrum bi^i 2 i 3 asEJ 



4?r (hhl 



(2Zi + 1)(2Z 2 + 1)(2Z 3 + 1) V° 



'3 
-1 



El 1 2 3 ) (ai 1 m 1 ai 2m2 a i 3 m 3 ) , (7) 
V mi rn,2 1^3 I 



mim2m3 
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Figure 1. The equilateral configuration of the reduced bispectra, \l 2 (I + l) 2 b;;i/(27i") 2 j X 10 16 . 
The weight is motivated by the behavior of the primary component, {Primary ^ g^jj [ mc 

shows the primary bispectrum from inflation with /jvi = 100 (see eq. m))- Dashed line shows the 
Sunyaev-Zel'dovich effect and weak lensing coupling in the Rayleigh-Jeans limit, while dotted line 
shows the extragalactic radio and infrared point sources. 



where the matrix is the Wigner-3j symbol. bi ± i 2 i a is found to be more useful to de- 
scribe the physical properties of the angular bispectrum than the conventional an- 
gular averaged bispectrum, B lMs = £ mim , m3 ( J, 2 ™ ) ( a hmi a hm 2 ai 3m3 ) . 

/ ' i 2 ''''3 / 

Note that Magueijdid defines a similar quantity, Bi x i 2 i 3 — bi^^/VAn. n 
We give the explicit form of the primary reduced bispectrum in the referencea. 
Figure 1 plots the equilateral case (h = 1% = h = reduced bispectrum. The 
acoustic oscillation in gTi{k) appears in bin. For the scale invariant <I>£-field power 
spectrum, 

b „y ^ 10 -17 /jvi x j-4 (8) 

apart from the oscillation and the Silk damping. 
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2 The CMB non-Gaussianity from low-redshift universe 



Secondary bispectra from the coupling between the Sunyaev-Zel'dovich (SZ) effect 
and the weak lensing effecttlHl3, and the extragalactic radio and infrared sourcesali3 
are found to be marginally detectable by MAP, while Planck will detect them 
significantly. 

We stive the explicit form of the SZ-lensing coupling bispectrum in the 
referenced. It is roughly 

jff-len, „ 10 -19 ju T p0 b gas X Z" 3 , (9) 

up to I ~ 1000. T P Q ~ 0.3 is the present day density weightecLmean temperature 
in units of 1 keV, while b gas ~ 6 is the linear gas pressure biasEJ. j v isJie spectral 
function of the SZ effect, and equals —2 in the Raylcigh Jeans regimeEl 

The CMB bispectrum from extragalactic point sources is easy to estimate, as 
long as they are assumed to be Poisson distributed. The reduced bispectrum is 
just constant in that case, i.e., &f S ;, ; = b ps = constant. Detailed arguments are 
summarized in the reference!. The results are 

b ps (90 GHz, < 2 Jy) - 2 x 10" 25 , (10) 
b ps (2l7 GHz,< 0.2 Jy) ~ 5 x 10" 28 . (11) 

Equations ([!(]) and ( |TT| ) correspond to MAP and Planck experiments, respectively, 
after subtracting > 5er sources from the map. 

Fi gure 1 plots ens in the Ray leigh— Jeans limit and b ps . Clearly, secondary 
sources dominate the primary signal as going to smaller angular scales, even if 
Jnl = 100 is taken. This is opposite to the CMB power spectrum that is dominated 
by either the primary CMB or the instrumental noise at all I in the COBE, MAP, 
and Planck experiments. Thus, we need to separate the primary bispectrum from 
secondary bispectra. Fortunately, the shapes of bispectra are fairly different, so 
that we can discriminate among them based on shapes. 



3 Can we detect /separate the CMB bispectra? 

Since the CMB bispectrum is very weak compared to the Gaussian noise from CMB 
itself, it is totally hopeless to measure &; 1 ; 2 ; 3 per mode. Instead, we fit the observed 
data by theoretical bispectra, and then constrain parameters. For this purpose, we 
minimize \ 2 given by 

x , (Ai)s (12) 



2<h<h<h 



where i denotes a component such as the primary, the SZ and lensing effects, extra- 
galactic sources, and so on. In case that the non-Gaussianity is small, the ccj 
variance of the bispectrum of i i is given by the six-point function of ai % 
of I l is dominated by either the primary CMB or the instrumental noise, while 
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Table 1. The signal-to- noise ratio for the detection of bispectra. /jvz, is the non-linear coupling 
parameter (Eq. (la)). T p o ~ 0.3 is the present day density weighted mean temperature in units 
of 1 kcV, while b gas ~ 6 is the linear gas pressure bias. j„ is the spectral function of the 
Sunyaev-Zel'dovich effect that has —2 in the Rayleigh-Jeans regime, b^ = b ps /10~ 25 , and 
627 = b pB /W~ 27 , where bP" is the reduced Poisson bispectrum. 
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Table 2. The signal degradation parameters, is defined so as di = 1 for perfectly separated 
component, while di > 1 for poorly separated one. 





primary 


SZ-lensing 


point sources 


COBE 


1.46 


3.89 


3.45 


MAP 


1.01 


1.16 


1.14 


Planck 


1.00 


1.00 


1.00 



the secondary sources are subdominant as long as COBE, MAP, and Planck exper- 
iments are considered. Taking dx 2 /dAi = 0, we obtain the Fisher matrix, 

^ n{>) 

bij = 2^ — ^2 ■ ( 13 ) 

2<h<l 2 <h 'i' 2 ' 3 

Since the covariance matrix of Ai is , we define the signal-to- noise ratio 

(S/N)i = 1 / y F~ 1 for a component i, and the signal degradation parameter 

di = FaF^ 1 due to a contamination from correlated components, di is defined 
so as di — 1 for perfectly separated component, while di > 1 for poorly separated 
one. Table 1 and 2 summarize (S/N)i and di, respectively, for COBE, MAP, and 
Planck experiments. In order to obtain (S/N) pr i mary > 1, /jvl > 600, 20, and 5 
are needed for COBE, MAP, and Planck experiments. As derived in section 1.1, 
the single field slow-roll inflation generally predicts Jnl ~ O(10~ 2 ), while the sec- 
ond order gravity yields Jnl ~ 0(1)- Therefore, either way, the detection is not 
expected for those experiments. 

(S/N) sz _i ens ~ 1 and (S/N) pa - 4 are expected for MAP, and (S/N) sz _ lens - 
20 and (S/N) ps - 30 for Planck. Thus, the leading contribution to the CMB 
bispectrum would be from extragalactic point sources. Then, the SZ- lensing con- 
tribution would be measured by Planck, yielding a constraint on T p0 b gas . This 
quantity measures the physical state of intergalactic gas pressure that is uncertain 
at present. 

4 Tentative comparison to the COBE 4 year bispectrum 

Ferreira, Magueijo, and GorskiEl claimed a significant detection of the equilateral 
bispectrum at I — 16 in the publicly available COBE 4 year data. Banday, Zaroubi, 
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and Gorskllj found that their detection would be ascribed to the systematics as- 
sociated with the eclipse of COBE satellite. This partly demonstrates that the 
bispectrum is sensitive to the non-linear effect even including the systematics in 
data. J—. 

MagueijoEJ tabulates a particular configuration of the COBE bispectrum with- 
out eclipse data. He tabulates 8 modes of 

Jf EE / bl ~ m+1 , (14) 

for I = 4 — 18. Figure 2 compares our primary bispectrum to his data. There is no 
significant detection of Jf per mode, and then /jvl is constrained as Jnl < 10 4 . It 
should be stressed that only 8 modes were used to constrain /jvl here. Since the 
total number of modes available in the full COBE bispectrum is 466 up to I = 20, 
the constraint on Jnl becomes much tighter than < 10 4 , say, la ~ 600 as suggested 
by our estimation (table 1). MAP and Planck will certainly put a very stringent 
constraint on /jvx, and thus test the inflationary scenario in non-linear order. 

5 Discussion 

What if the primary bispectrum is detected? In other words, what if /nl 3> 1 
is measured? One might suspect the systematics, as already demonstrated in the 
COBE data, or foreground secondary sources. Since we fit the theory to the data, 
the goodness- of- fit value of x 2 a t the best-fit /jvl would tell us whether detected 
signal is really primary or not. If it passes the % 2 test, then the most straightforward 
consequence is that — the single field slow-roll inflation is not right. 
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